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Abstract.  Solar  cells  have  been  made  from  hydrogenated  amorphous  silicon  (a-Si:H)  films.  The  films  including  Si  nanopar¬ 
ticles  by  3%  in  the  volume  fraction  have  been  found  to  possess  good  properties.  In  order  to  find  the  method  to  control  this 
volume  fraction,  the  motion  of  Si  nanoparticles  in  the  rarefied  flow  of  H,  and  SiH4  is  examined  for  a  CVD  (chemical  vapor 
deposition)  reactor  by  use  of  the  DSMC  method. 


INTRODUCTION 

If  every  house  in  the  world  had  a  wide  solar  cell  on  its  roof,  a  large  reduction  of  CO-,  around  the  globe  would  surely  be 
realized.  Solar  cells  are  made  from  hydrogenated  amorphous  silicon  (a-Si:H)  films  or  microcrystalline  films.  Recently, 
a  better  alternative  is  proposed  by  the  group  of  Kyushu  University  [1].  It  is  the  a-Si:H  films  into  which  nanocrystalline 
particles  of  silicon  are  dispersed.  This  film  has  the  light  absorption  coefficient  as  large  as  the  usual  a-Si:H.  Moreover, 
it  shows  a  small  light-induced  degradation  and  a  good  carrier  transport  properties  as  in  the  case  of  the  microcrystalline 
films. 

In  order  to  realize  the  optimum  property  for  such  nanoparticle-dispersed  films,  the  transport  of  Si  nanoparticle  to 
the  surface  of  deposited  a-Si:H  films  must  be  controlled  in  such  a  way  that  the  volume  fraction  of  Si  nanoparticles 
dispersed  into  the  film  is  3%.  To  find  a  method  to  control  the  volume  fraction,  the  motion  of  the  nanoparticles  in  rarefied 
gas  flows  of  H,  and  SiH4  is  examined  for  a  CVD  (chemical  vapor  deposition)  reactor  using  the  DSMC  simulation. 


MODELING  OF  PHENOMENA 

Figure  1  shows  the  CVD  reactor  treated  here.  There  are  two  electrodes  in  the  cylindrical  chamber  with  diameter  d3  and 
height  /;,.  The  diameter  and  height  are  dx  and  for  the  upper  electrode,  and  d2  and  /;  ,  —  h2  for  the  lower  electrode. 
The  upper  electrode  is  powered  by  a  radio  frequency  source  and  the  lower  is  grounded.  Source  gases  H,  and  SiH4  are 
fed  uniformly  from  the  shower  head  set  on  the  lower  surface  of  the  upper  electrode.  Nanoparticles  of  silicon,  (Si),,,  are 
generated  by  the  electric  discharge  between  the  two  electrodes;  SiH4  molecules  are  dissociated  by  electron  impact  and 
released  Si  atoms  grow  into  a  nanoparticle  by  repetition  of  collisions  with  other  Si  atoms.  Gases  and  nanoparticles  are 
pumped  out  from  the  annular  exist  at  z  =  h,. 

The  mechanism  of  the  generation  and  growth  of  a  nanoparticle  is  complicated.  Replacing  the  mechanism  by  a  very 
simple  model,  we  examined  the  flows  of  H,,  SiH4,  and  (Si)„  in  the  reactor  of  Fig.  1.  There  is  a  great  difference  of 
molecular  mass  between  the  carrier  gas  and  nanoparticle.  Our  concern  is  the  velocity  slip  among  H,,  SiH4,  and  (Si)„. 


OUTLINE  OF  SIMULATION 

Gas  is  a  mixture  of  Ho,  SiH4,  and  (Si)„.  Rarefied  flows  of  the  gas  mixture  are  calculated  using  the  DSMC.  A  three- 
dimensional  cell  network  is  introduced  for  a  general  use.  The  cell  sizes  are  A r  and  Az  in  the  r-  and  z-  direction, 
respectively.  They  are  constant.  The  cell  boundary  (grid  point)  is  given  by  r(  =  iAr  and  z  ■  =  jAz.  The  region  between 
r_  j  and  r.  has  the  area  (2 i—  l)n(Ar)2.  The  region  is  divided  into  K[=  2(2 i—  1)]  subregions  in  the  8-  direction.  The 
cell  boundary  is  8k  =  kAO,  where  k  =  1,2,  -  ■  -  ,K  and  AO  =2 tz/K.  The  cell  volume  Vc  is  Tt(Ar)2Az/2  for  any  cell.  Free 
motion  and  collision  of  simulated  molecules  are  treated  separately  in  the  DSMC. 
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FREE  MOTION 


Let  1,  2,  and  3  denote  the  species  H2,  SiH4,  and  (Si)„,  respectively.  Since  there  is  a  large  difference  among  the  number 
densities  of  three  species,  we  introduce  weight  Ws  for  a  simulated  particle  of  species  s,  where  s  =  1,  2,  or  3.  Let  Ql  and 
Q7  be  the  mass  flow  rates  of  H-,  and  SiH4  in  units  of  seem.  These  gases  are  fed  from  the  shower  head.  Since  H-,  is  a 
carrier  gas,  usually  Q1  Q2-  We  choose  W2  from  W,/W]  =  Q-,/Qx.  The  weight  Wl  is  determined  from  Wt  =n0Vc/NQ, 
where  nQ  is  the  number  density  of  H2  and  N0  is  the  number  of  simulated  particles  in  a  cell  at  time  zero.  The  number  of 
incoming  molecules  in  At  is 

p  _  nSTpQiAt 
'  60  x  106W,. 

where  nSTp  is  the  standard  number  density  at  273  K  and  1  atm.  Our  choice  of  Wt  yields  F]  =  F2. 

The  nanoparticle  (Si)„  is  assumed  to  be  produced  by  the  reaction. 


«SiH4  — »  (Si),,  +  2nH2  (1) 

Let  mQ  and  m3  be  the  mass  of  Si  and  (Si)„.  We  see  that  the  number  of  SiH4  molecules  in  producing  one  (Si)„  is  m 
The  mass  flow  rate  of  SiH4  is  Q*,.  A  part  of  Q-,  is  consumed  for  reaction  (1)  and  the  rest  of  Q2  is  pumped  out  from  the 
exit.  Let  T]Q0,  where  0  <  rj  <  1,  be  consumed  in  the  reaction.  Therefore, 


os  =  ^~ 

m3/mo 


(2) 


is  the  rate  of  production  of  (Si);i.  It  is  assumed  that  (Si)„  is  produced  spatially  uniformly  in  the  region  between  the 
electrodes.  We  choose  W3  as  W3/W2  =  mQ/m3.  If  77  =  1  in  eq.(2),  this  means  that  W2/W2  =  Q2/Qt- 

The  shower  head  has  many  tiny  holes.  The  flow  of  the  mixture  of  H  ,  and  SiH4  is  regarded  as  sonic  at  the  exit  of  the 
hole.  That  is,  the  flow  velocity  V*  and  temperature  T*  are  given  by  [2] 


V* 


* 


s/yRT7 

2 

r+ir° 


822 
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where,  y  is  the  specific  heat  ratio  for  the  mixture,  T0  is  the  stagnation  temperature  that  is  close  to  the  room  temperature, 
R  =  k/m,  k  is  the  Boltzmann  constant,  and  m  is  the  molecular  mass  of  the  mixture.  The  velocities  of  incoming 
molecules  are  sampled  from  the  Maxwellian  distribution  with  the  flow  velocity  V* .  The  gas  pressure  between  the 
two  electrodes  is  set  to  a  given  value  by  adjusting  the  transmission  probability  of  a  particle  at  the  annular  exit.  The 
diffuse  reflection  model  is  used  on  any  solid  wall. 


MOLECULAR  COLLISION 


As  before,  let  1,  2,  and  3  denote  H2,  SiH4,  and  (Si)„.  There  are  six  types  of  collision.  Like  collisions  are  1-1,  2-2, 
and  3-3,  and  unlike  collisions  are  1-2,  1-3,  and  2-3.  The  rigid  sphere  model  is  used  for  three  species.  Let  <ys  be  the 
molecular  diameter.  We  set  ctj  =  0.2745  nm  and  a-,  =  0.4187  nm,  where  the  diameter  of  CH4  is  used  instead  of  the 
diameter  of  SiH4  [3],  The  diameter  dn  of  nanoparticles  is  0.5  —  5  nm  [1],  Here  dn  is  the  input  data.  The  mass  m3  of  a 
nanoparticle  is  chosen  as 


where  p(=  2.33  x  103kg/m3)  is  the  density  of  solid  silicon  in  bulk.  We  have  m3/m2  =  22.90  for  dn  =  1  nm.  The 
collisions  in  each  cell  are  treated  using  the  maximum  collision  number  method  [4],  Let  /vj(,nax)  be  the  maximum 
collision  number  between  species  s  and  r.  It  is 


jy(max)  _  q  (Ns  8sr)Nr 


Vc 


(3) 


where 

Csr  =  — max(Ws ,  Wr)g^x]  (4) 

1  H-  osr 

where  8sr  is  the  Kronecher  delta,  Ns  is  the  number  of  simulated  molecules  of  species  s  in  a  cell,  the  diameter  <7sr  is 
defined  by  n(ds  +  dr)2 /4,  and  is  the  maximum  value  of  the  relative  speed  between  two  molecules  in  species  s 

and  r.  It  is  estimated  as  =  2.5  y/kTref/ psr,  where  jlsr  =  msmr/ ( ms  +  mr),  and  ^ref  is  a  reference  temperature.  The 

procedure  to  calculate  s-r  collisions  in  a  cell  is  as  follows  [4],  Repeat  the  following  steps  (i) — (ii)  /vj™x)  times. 


(i)  Sample  two  particles  randomly,  one  from  species  ,v  and  another  from  species  r.  Let  (cs,  c,-)  be  the  velocities  of 
this  pair.  Calculate  gsr  =  |cs  -  cr|.  If  U  <  gsr/g[™x\  go  to  (ii),  and  if  not,  go  to  the  next  do-loop  of  A^max) .  Here 
U  is  a  uniform  random  number  of  0  <  U  <  1 . 


(ii)  Replace  cs  by  c'.  with  probability  Wrj max(Ws,  Wr)  and  cr  by  c'r  with  probability  Ws/ max(Wy,  W,-),  where  (c's,  cj.) 
are  the  post-collision  velocities  given  by 

c',  =  W  +  —gsrR 

s  M° 


c  =  W- 


M 


i'srK 


where  M  =  ms  +  mr,  W  =  (mscs  +  mrcr)  / M ,  and  R  is  a  unit  vector  with  a  random  direction. 


RESULT  AND  DISCUSSION 

Computational  conditions  are  as  follows.  The  gas  pressure  p  of  H2,  which  is  close  to  the  total  pressure,  is  10  —  20 
Pa,  the  mass  flow  rate  Qx  is  500  —  1000  seem,  Q2  =  Ql/500,  and  the  diameter  dn  of  nanoparticle  is  0.5  —  2  nm.  The 
geometry  of  the  CVD  reactor  is  dx  =  d0  =  80  mm,  d3  =  100  mm,  h  ]  =30  mm,  h2  =  50  mm,  and  /;  ,  =  90  mm. 

First  we  show  the  flow  fields  of  H2,  SiH4,  and  (Si)„  for  Qj  =  500  seem,  Q0  =  1  seem,  p  =  20  Pa,  and  dn  =  1  nm. 
Figures  2,  3,  and  4  show  the  flow  velocities  of  three  species.  The  flow  of  SiH4  is  similar  to  that  of  H2.  Two  vortices 
exist  near  the  corners  of  the  upper  and  lower  electrodes.  Nanoparticles  flow  slowly.  The  flow  of  (Si)„  near  the  edge  of 
the  electrode  is  similar  of  that  of  the  carrier  gas  H2  but  the  lower  vortex  does  not  appear  for  (Si)„.  The  flow  of  (Si)„ 
near  the  exit  is  very  slow,  almost  retarded.  Figures  5,  6,  and  7  show  the  number  densities  of  H2,  SiH4,  and  (Si)„.  As  for 
H2,  the  density  decreases  in  the  downstream  region.  But  for  SiH4  and  (Si)„  the  opposite  is  true.  Figures  8,  9,  and  10 
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FIGURE  4.  Velocity  of  nanoparticle  (dn  =  1.0  nm)  FIGURE  5.  Density  of  H2 


show  the  effect  of  the  size  dn  on  the  flow  velocity  of  nanoparticles.  The  lower  vortex  found  for  dn  =  0.5  nm  disappears 
for  dn  =  2  nm.  Figures  11,  12,  and  13  show  the  number  densities  for  dn  =  0.5,  1.0,  and  2.0  nm.  For  all  three  cases  the 
density  is  highest  near  the  exit.  The  density  gradient  along  the  streamline  from  the  shower  head  to  the  exit  increases 
with  the  size  of  nanoparticle.  We  examined  the  effect  of  pressure  on  the  velocity  and  density  fields  of  nanoparticles  by 
comparing  the  data  for  p  =  10  Pa  and  20  Pa.  We  have  seen  that  the  effect  is  small.  Also,  we  examined  the  effect  of  the 
mass  flow  rate  on  the  flow  field  of  nanoparticles,  by  comparing  the  data  for  Qx  =  500  and  1000  seem.  The  pressure  is 
20  Pa  and  the  size  dn  is  1  nm.  Although  the  velocity  is  higher  for  1000  seem,  the  effect  on  the  density  field  is  small.  In 
order  to  fabricate  nanoparticle-dispersed  films  with  the  best  quality,  control  of  the  flux  of  nanoparticles  onto  the  lower 
electrode  is  important.  Figure  14  shows  the  flux  for  Q  t  =  500  seem,  Q2  =  1  seem,  and  p  =  20  Pa.  We  see  that  the  flux  is 
uniform  on  the  electrode  in  particular  for  dn  =  0.5  and  1.0  nm.  This  means  the  correctness  of  the  choice  of  the  present 
operating  condition.  The  flux  is  approximately  proportional  to  d~ 4.  The  effect  of  Ql  on  the  flux  is  examined  for  p  =  20 
Pa  and  d„  =  1  nm.  The  effect  is  small;  The  flux  for  Q1  =  1000  seem  is  larger  by  10%  than  the  flux  for  Ql  =  500  seem. 
Figure  15  shows  the  effect  of  pressure  p  on  the  flux.  When  p  is  halved  from  20  Pa  to  10  Pa,  the  flux  is  also  halved; 
The  flux  is  proportional  to  p. 
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FIGURE  6.  Density  of  SiH4  FIGURE  7.  Density  of  nanoparticle  (dn  =  1.0  nm) 
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FIGURE  12.  Density  of  nanoparticle  (dn  =  1.0  nm) 


FIGURE  13.  Density  of  nanoparticle  (dn  =  2.0  nm) 
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FIGURE  14.  Flux  of  nanoparticles 


FIGURE  15.  Flux  of  nanoparticles 


CONCLUSIONS 

Flow  field  of  silicon  nanoparticles  in  the  rarefied  flow  of  H  ,  and  SiH4  is  examined  using  the  DSMC.  The  pressure  p 
is  10  —  20  Pa,  the  mass  flow  rate  Ql  of  H2  is  500  —  1000  seem,  and  the  size  dn  of  nanoparticle  is  0.5  —  2.0  nm.  We 
conclude  as  follows. 

1.  The  flux  of  nanoparticles  onto  the  wafer  is  nearly  proportional  to  d~4. 

2.  The  flux  is  almost  independent  of  Qx . 

3.  The  flux  is  nearly  proportional  to  p. 
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